Supercurrent through grain boundaries in the presence of strong correlations 
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Strong correlations are known to severely reduce the mobility of charge carriers near half-filling 
and thus have an important influence on the current carrying properties of grain boundaries in the 
high-Tc cuprates. In this work we present an extension of the Gutzwiller projection approach to 
treat electronic correlations below as well as above half-filling consistently. We apply this method to 
investigate the critical current through grain boundaries with a wide range of misalignment angles 
for electron- and hole-doped systems. For the latter excellent agreement with experimental data is 
found. We further provide a detailed comparison to an analogous weak-coupling evaluation. 

PACS numbers: 74.81.-g, 74.25. Sv, 74.50.+r, 74.45. +c 



Introduction. — The manufacturing of high- 
temperature superconducting cables as well as the 
characterization of well-defined interfaces and contacts 
for superconducting microwave electronics have been 
optimized over the last 20 years so that they are now 
found in applications ranging from short range power 
supply to medical instrumentation. This was made 
possible by intensive experimental and theoretical 
research, which has led to a better understanding of 
the bulk as well as the interface properties of high-T c 
materials One of the most complex and puzzling 

problems in this respect is the strong reduction of the 
critical current function of the grain boundary 

(GB) angle, which has been extensively studied on 
samples with artificially fabricated, well defined GBs 
Although great progress has been made to reduce the 
influence of large angle GBs on the total current [l(| and 
to locally improve their current carrying properties 
a full theoretical understanding of the grain boundary 
problem is hampered by the complexity of the physics 
involved. In a recent study the experimentally observed 
decay of j c with increasing misalignment angle was 
reproduced in a microscopic model in which charge 
inhomogeneities at the grain boundary were identified 
as the main source for the suppression of j c 12]. But 
still the overall magnitude of j c could not be correctly 
determined and it was speculated that strong electronic 
correlations present in the high-T c materials might be 
responsible for this discrepancy jl3|. 

In this article we present an analysis of the dependence 
of the supercurrent on the grain boundary angle in the 
presence of strong Coulomb interactions. In order to de- 
scribe the generic lattice defects for a certain GB angle, 
a theoretically reconstructed sample of a GB should in- 
clude at least several hundred lattice sites. Due to numer- 
ical limitations, this can only be modeled by an effective 
one-particle description. To incorporate static correla- 
tions nevertheless, we employ a Gutzwiller approxima- 
tion [Tij that should already capture an important part 
of the GB physics. The Gutzwiller approach has had 
considerable success for homogeneous models of cuprate 
high-r c superconductors [l5| . Only in recent years was it 
rigorously extended to inhomogeneous systems [HI, E3] , 
but is now an established method to treat important as- 
pects of the interplay between impurities and supercon- 



ductivity [18|] . 

The model system of a GB is, in a certain respect, a 
particularly inhomogeneous system that deals with very 
strong impurities that to our knowledge have not yet 
been treated within the Gutzwiller approach. The en- 
ergy scale of the system is set by the hopping parame- 
ter £t>uik in the homogeneous bulk regions far away from 
the GB. In cuprates we typically expect a Coulomb re- 
pulsion of U ~ 20 tbuik and a spin exchange interaction 



Jij ~ ta, where t 



ij, wiicic is the hopping matrix element for 



nearest and next-nearest neighbor sites. It has been 
shown that charge barriers at the GB lead to local impu- 
rity potentials of the order of £j ~ — 20ibuik [13] or more, 
therefore close to the GB, the estimate £j ~ —U applies. 
Although in the bulk regions of a hole-doped system, the 
local charge densities nt are below half filling, the pos- 
itive charges in the vicinity of the GB can induce sites 
with m well above half filling, even in the presence of 
strong Coulomb repulsion. As an example, in Fig. [He), 
(d) and (e) local charge densities, impurity potentials and 
hopping amplitudes are shown as a function of the dis- 
tance from a GB with misalignment angle a — 44° . Such 
systems with strongly positively charged regions cannot 
be described by the standard Gutzwiller approximation. 
In the first part of this article we present an extended 
Gutzwiller formalism appropriate to describe such sys- 
tems. In the second part we analyze our results for the 
critical current calculated with the new method and show 
that they are in excellent agreement with experiments. 

Model and methodology. — We use the Gutzwiller 
projected t-J model to describe the GB system. Within 
the Hartree-Fock decoupling scheme it is given by the 
Hamiltonian 



ijs 



h.c. 



with Ay = (|4 + i) JyAy, 
Xij = (j9ij ~ j)JijXiji 



(1) 

(2) 
(3) 



where Ay = \{{c liCn ) + {c u c lt }), Xij = |((4t c it) + 

( c li c Ji))' ™i = J2s c ls c is> and Hi = V - £i- The matrix 
for the tunneling amplitudes Uj has non-zero entries for 



2 



nearest and next-nearest neighbor hopping. The tunnel- 
ing amplitudes and the local potentials e$ are calcu- 
lated using the molecular dynamics algorithm of Ref. [ID 
and reflect the charge barriers and defects of the lat- 
tice at the GB. The spin-interaction coefficients Jy are 
obtained by rescaling the nearest-neighbor entries of fy 
with a constant factor. 

The Gutzwiller factors in ([T]) are defined as follows: 



g\g\ and g 



J ~~ 9i9'\i where 



9t = 



2\l-rn\ 

II-THI + I' 

2 



1-Tli +1 



(4a) 
(4b) 



For local densities rii < 1 these expressions coincide with 
the commonly employed definitions of Ref. By in- 
troducing the absolute values in their definitions we ob- 
tain a renormalization that is symmetric with respect to 
half filling rii = 1 . This is equivalent to a local particle- 
hole transformation for sites where nt > 1. Due to the 
particle-hole invariance of the Hubbard model from which 
the t—J model is derived, this treatment is obviously cor- 
rect for a homogeneous system with rii = n > 1. 

The rigorous derivation for an inhomogeneous sys- 
tem is based on the definition of a projection operator 
V = Yli 'Pi that locally projects either doubly occupied 
or empty sites out of a BCS type wave function \ip)o, 
depending on the local density: 



yf*(l - A) if «i < 1 
yfil-Ei) ifn,>l 



(5) 



Here, yi is a fugacity factor, while Di — h^n^ and 
Ei = (1 — hi^)(l — hn) denote the projection on dou- 
bly occupied and empty sites, respectively. Note that 
the selection of different projectors depending on the lo- 
cal density is a natural extension of the operator y?* : this 
operator assigns different weights to contributions from 
empty and singly occupied sites and thus ensures a pro- 
jected state \tp) that can be transparently compared with 
the pre-projected state |^>o) due to its common set of lo- 
cal densities [T(|. A comparison of expectation values 
evaluated with the wave functions \tpo) and \tp) = V\tjM), 
completely analogous to the one presented in Ref. Ila . 
yields the results of Eq. ((4]) . 

It remains to be justified that the implementation of a 
symmetric renormalization around half filling, or equiv- 
alcntly, projecting out empty and doubly occupied sites 
in the free wave function, is physically meaningful. The 
argumentation here is very similar to the original idea 
of Gutzwiller, to project doublons out of the wave func- 
tion due to their energy cost. Our extension is obviously 
justified for bulk regions: In an environment where [ii is 
large enough to yield n» > 1, a local particle- hole trans- 
formation leads to the removal of empty sites from the 
wave function. Also at the boundary of such a region 
the definition of Eq. ([5|) is consistent since on a site with 
Hi ~ U, the gain of potential energy cancels the cost of 
double occupancy, independent of the situation on neigh- 
boring sites. Thus the energy cost to create a holon at 
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FIG. 1. The lattice structure for a (520) GB with misalign- 
ment angle a( 52 o) = 44° in the a;,y-plane (a). Orange lines 
mark the bonds defining the "channels" through the GB that 
explicitly appear in the sum of Eq. |7|l. Results of a self- 
consistent calculation for the (520) GB at a bulk filling of 
n — 0.8 without (i) and with (ii) Gutzwiller factors include 
the d-wave pairing amplitude Aij (b) and the local density 
rii (c). In addition the local impurity potentials (d) and 
the averages over all nearest neighbor (t) and next-nearest 
neighbor (t') hopping parameters of single sites i (e) are 
shown. Energies are measured in units of tbuik = 1. The 
spin-interaction is calculated from the nearest-neighbor en- 
tries of tij for (i) as Jij = 2.38tij and for (ii) as Jij = 0.9tij 
so that Abuik = 0.275£buik has the same value for both calcu- 
lations. The precise values of the impurity potentials Ei are 
obtained from the charge fluctuations assuming a screening 
length I = 2 A. 



site i is of the size 
out. 



of U and consequently it is projected 



Note that this procedure is possible in the Gutzwiller 
approximation only due to its fundamental assumption 
that inter-site correlations can be neglected. However, in 
the framework of the Gutzwiller approach it is a consis- 
tent improvement of the wave function modeling in the 
spirit of assigning different weights to different sites with 
the fugacity factor y"*. It is a heuristic, very detailed and 
sophisticated method of minimizing the energy of a sys- 
tem described in a phenomenological mean-field picture. 
Below we will show that this picture is sufficiently realis- 
tic to accurately describe the physics of superconducting 
charge transport through a GB. 

The renormalized model Hamiltonian of Eq. (JlJ is eval- 
uated within the Bogoliubov-de Gennes (BdG) frame- 
work. We enforce constant phases with a difference 
<fi = tt/2 on the pairing amplitude Ay on either end 
of the sample, in the center of which the GB is situated. 
The corresponding phase gradient induces finite currents 
given by 



Jij 



.))• 



(6) 



3 



The normalized total current through the GB is 
sin(a/2) 



.7c 



E 



H 



(7) 



Xi<0,Xj>0 



where Xi and Xj are the x-coordinates of sites on either 
side of the GB, a is the misalignment angle and M y the 
second Miller index of the (M x M y M z ) GB. The factor 
sm(a/2)/ My normalizes the current to a unit length in 
y-direction. To be able to compare with experimental 
results we employ the following relation [12[ between the 
current j^ bs in absolute units A/cm 2 and the current in 
units [ibuik] = eV 



■abs 



Jc 



A cx p N UG e 



v thco 



ac 



(8) 



where A^uc = 2 is the number of Cu02 planes per 
unit cell, a the misalignment angle and a = 3.82A and 
c = 11. 7 A the lattice spacing in x- and z-direction, re- 
spectively. The values A ox p - 0.025eV and A thco = 
Abuik = 0.275eV are valid for the bulk regions. In order 
to compensate the smaller magnitude of the experimental 
quantities we employ the factor A exp /A thco ~ 0.09 where 
we make use of the fact that j c scales approximately lin- 
early w.r.t. A if A is small. We cannot employ the ex- 
perimental value A oxp ~ 0.025eV in our calculations as 
for gap values of this order of magnitude, the iterative 
solution of the BdG equations is non-convergent. 

Angle dependence of the super current. — In Fig. [2(a) 
and (b) we show on a logarithmic scale the dependence of 
the zero-temperature supercurrent on the GB misalign- 
ment angle for a hole-doped system. We compare two 
types of calculations with experimental data taken from 
Ref. U (i) a standard evaluation of the GB-Hamiltonian 
within the BdG framework (i.e. g\^ = gf^ = 1) and (ii) 
an evaluation with the local Gutzwiller factors as defined 
in Eq. (fij). In order to provide a meaningful comparison 
of (i) and (ii) we employ the same hopping matrices Uj 
and impurity potentials £j for all calculations, but glob- 
ally scale Jij to obtain identical values for the d-wave 
pairing amplitude in the bulk: Abuik = 0.275 ibuik- This 
approach is appropriate as j c tx Abuik and it is in com- 
plete analogy to Ref. [H. Note that also at the GB the 
pairing amplitude A y - shown in Fig. [lja) does not differ 
qualitatively for the evaluation without (i) and with (ii) 
Gutzwiller factors. 

As can be seen from Fig. [2] the supercurrent decays in 
both cases (i) and (ii) exponentially with the GB angle. 
We also find that the current in the correlated system 
(ii) is almost one order of magnitude smaller than in the 
system without correlations (i). As the screening length 
I in the Cu02 planes is only approximately known we 
performed two calculations assuming different values of 
I, with I = 2k for Fig. [2£a) and I = 1.2k for Fig. 02b). 
For the choice of I = 2k we obtain excellent agreement 
with the experimentally determined critical currents over 
a wide range of GB angles. 

In Fig. [21(c) we show the angle dependence of the cur- 
rent for an electron-doped system. Comparing panels 
(b) and (c) of Fig. [2l we find that the current is almost 
particle-hole symmetric. Correlation effects are known 
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FIG. 2. Angular dependence of the critical current for GB 
samples with electron density n = 0.8 (a,b) and n = 1.2 
(c). We plot results for 11 samples for 6 different angles in 
between 15° and 45° in order to be able to average. This av- 
erage is calculated as a linear fit to the logarithm of the data 
points and depicted as colored lines. The impurity potentials 
have been calculated from the charge fluctuations assuming a 
screening length of I = 2A (a) and I — 1.2A (b,c). The cur- 
rent is given in units of tbuik- Depicted are calculations using 
the BdG Hamiltonian of Eq. (JXJ) where for (i) all Gutzwiller 
factors are set to 1 and for (ii) all Gutzwiller factors are taken 
as defined in Eq. Q. The theoretical values in (a) and (b) 
are compared to experimental data (iii) for optimally doped 
YBCO thin films on SrTi0 3 taken from Ref. 1. The spin- 
interaction is calculated from the nearest-neighbor entries of 
tij for (i) as Jij = 2.387^- and for (ii) as Jij — 0.9tij so that 
Abuik = 0.2757buik has the same value for both calculations. 



to be strongest at half filling where a Mott insulator is 
obtained in the Gutzwiller approximation. The domi- 
nant reduction of the current due to correlations should 
therefore be caused by areas in the system where n» ~ 1 . 
These areas are equally present in the hole- and electron- 
doped system in the transition region to the bulk. In 
the hole-doped case such a region is already implied by 
the fact that, directly at the GB, the system is always 
well above half filling. In the electron-doped case it is 
caused by large density fluctuations. These appear due 
to the charge fluctuations with the same magnitude in 
the electron- and the hole-doped case. For a hole-doped 
example, Fig. [ljc) shows values for the density with 
0.5 < ni < 2. Obviously, these arguments imply that 
the correlation effects depend on the charge fluctuations. 
Furthermore, the fact, that the narrow transition regions 
have the same extension for all angles, explains that the 
exponential behavior w.r.t. the GB angle appears in both 
calculations (i) and (ii): only the break-up of pairs by i 
mpurities, present in (i) and (ii), causes the exponential 
behavior while the correlation induced reduction of the 
supercurrent is constant for all angles. 

Fig. [3] displays the microscopic (local) values of the 
currents jij, the Gutzwiller factors g\, and the local po- 



4 



10 
L(f 
10"' 
s, 10 

* io" : 

10" ! 

io- 




X — X (i) without correl. 
H — h (ii) with correl. 

H 1 1 h 




2 4 6 8 10 12 14 16 
channel number 



Lr 
0.8- 



0.4- 

0.2- 



- i--I2A | 




(e)_ 


— :\) 













-50 



-25 



„ 25 
X[A] 



50 



1 2 3 4 5 6 7 
channel number 



FIG. 3. Cuts parallel to a (520) GB with angle a (520 ) ~ 44° 
in y-direction at x = — 12A (a,c) and x = (b,d,f). In the 
first row (a,b) we show the values for the currents jij through 
all channels of a cut. The channels of the x — cut are 
depicted in Fig. OJa) as orange lines. Negative currents are 
marked with a circle around the data point. In the second row 
(c,d) we show the Gutzwiller factors g\j and in (f) the local 
impurity potentials Si of the two sites that are connected by 
a channel. In panel (e) the average g l of gjj along cuts for 
fixed values of x is displayed. 



tentials e$ for cuts parallel to the GB. Panels (a,c) show 
these values at the edge of the GB (at x — — 12A) where 
potential fluctuations are minor but where the average g l 
of gjj along the cut is minimal, indicating the strongest 
suppression of hopping by Coulomb repulsion (see panel 
(e)). Panels (b,d,f) refer to a cut along the center of 
the GB (at x = 0), where the potential fluctuations are 
maximal but the Gutzwiller factors gjj are close to 1 
due to the strong effective electron doping induced by 
the local impurity potentials. The flat distribution of g\j 
in panel (d) leads to the homogeneous reduction for all 
channels in panel (b) while in panel (c) g\j varies signif- 
icantly which results in a non-homogeneous reduction of 
jij in panel (a). Current conservation guarantees that the 
total reduction of the supercurrent in (b) equals the to- 
tal reduction in (a). Both panels (a,b) show that charge 
transport is carried almost entirely by very few chan- 
nels. At the GB center these channels are exposed only 
to slightly varying impurity potentials in contrast to the 
ones that carry almost no current (panel (f)). Even at 
the edge of the grain boundary, where almost no poten- 
tial fluctuations are present, back currents continue to 
exist (panel (a)). They are induced by the impurities 
at the center of the GB and degrade slowly on account 
of the small phase gradient away from the GB. Finally 
we point out that fluctuations of the critical current are 
strong in experiment and that they are better reproduced 
by the correlated than by the non-correlated modeling 
(see Fig. |2Ja,b)). Random Coulomb blocking of some of 
the few channels that carry almost all current produces 
these correlation induced fluctuations. 

Conclusion. — We expect that the extension of the 
Gutzwiller projection approach, which we introduced to 



treat strongly inhomogeneous systems, is not only vi- 
able for the characterization of the supercurrent through 
grain boundaries. It should also be instrumental for other 
systems with intrinsic impurities or artificial structures 
(heterostructures) . We emphasize that we achieved con- 
vincing agreement with experiment for the supercurrent 
through GBs. The reduction of the current was shown to 
be due to the intertwining effects of large charge fluctua- 
tions and strong correlations. It can thus be understood 
that experimental doping of the GB in order to moder- 
ate these fluctuations significantly improves the current 
carrying properties, which has already been observed for 
Ca-doped interfaces [HI H3 • 
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